Random Access Broadcast: Stability and 
Throughput Analysis 

Brooke Shrader, Student Member, IEEE and Anthony Ephremides, Fellow, IEEE 



Abstract 

A wireless network in which packets are broadcast to a group of receivers through use of a random 
access protocol is considered in this work. The relation to previous work on networks of interacting 
queues is discussed and subsequently, the stability and throughput regions of the system are analyzed 
and presented. A simple network of two source nodes and two destination nodes is considered first. The 
broadcast service process is analyzed assuming a channel that allows for packet capture and multipacket 
reception. In this small network, the stability and throughput regions are observed to coincide. The 
same problem for a network with N sources and M destinations is considered next. The channel model 
is simplified in that multipacket reception is no longer permitted. Bounds on the stability region are 
developed using the concept of stability rank and the throughput region of the system is compared to the 
bounds. Our results show that as the number of destination nodes increases, the stability and throughput 
regions diminish. Additionally, a previous conjecture that the stability and throughput regions coincide 
for a network of arbitrarily many sources is supported for a broadcast scenario by the results presented 
in this work. 
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Random Access Broadcast: Stability and 
Throughput Analysis 

I. Introduction 

The stability and throughput of finite-user random access systems for unicast transmission have been 
studied extensively. The throughput analysis is included in some of Abramson's early work on the topic 
[1], while the stability problem was first introduced by Tsybakov and Mikhailov [2]. The finite-user 
stability problem proves to be much more difficult than the throughput problem, and as such, the history 
of the stability problem is more rich and interesting. In [2], sufficient conditions for ergodicity were found 
using the transition probabilities of the Markov chain corresponding to queue lengths. Later work [3] 
described the use of stochastic dominance as a means of characterizing the stability region. Additionally, 
stability conditions based on the joint queue statistics were provided in [4]. Despite these and other 
attempts, the stability region for arbitrarily many sources N remains unsolved. In [5] the concept of 
stability rank was introduced and provided the tightest known bounds to the stability region for N 
sources. The exact stability region for N sources and an arrival process which is correlated among the 
sources was obtained in [6]. 

Most works, including all of those mentioned above, study random access under the collision channel 
model, in which transmission by more than one source results in failed reception of all packets. Some 
recent works have incorporated the probabilistic nature of reception and the possibility of multipacket 
reception (MPR) into the channel model. The stability of infinite-user random access with MPR was first 
examined in [7]. More recently, the finite-user problem was examined in [8] and it was shown that the 
possibility of MPR results in an increase in the stable throughput of the system. The benefit to stability 
was so dramatic that for channels with sufficiently strong reception capabilities, random access was 
shown to outperform time division multiple access (TDMA) schemes. This result provides motivation for 
a renewed interest in random access. 

In this correspondence we introduce multiple destination nodes and broadcast transmission into the 
network and study the resulting stability and throughput performance of random access. The introduction 
of multiple destinations is a necessary first step in understanding the behavior of ad hoc and multihop 
networks, where random access presents an advantage over TDMA due to its distributed nature. In 
particular, we analyze the performance of a random access broadcast system in which a source node 
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Fig. 1. The broadcast scenario we study in this work: N source nodes transmit to M destination nodes. 



sends a common packet to all destination nodes. Broadcast transmission is useful for control of the 
network (ie, initialization, route discovery, timing synchronization) and for a number of applications. 

II. Model and formulation 

Consider a system S consisting of N source nodes, s%, S2, ■ ■ ■ , s/y, and M destination nodes S~ 2 \ . . . , 
Packets arrive to s n according to a Bernoulli process with rate X n , n = 1,2, ...,N packets per slot. 
The arrival process is independent from source to source and independent, identically distributed (i.i.d) 
over slots. Packets that are not immediately transmitted are stored in an infinite buffer maintained at each 
source. All source nodes compete in a random fashion for access to the channel in order to transmit 
a packet of information to all destination nodes. When source n has a packet to transmit, it does so 
with probability p n in the first available slot. This scenario is depicted in FigfT] Each packet is intended 
for all M destinations. We assume that instantaneous and error-free acknowledgements (ACKs) are sent 
from the destinations and that each source-destination pair has a dedicated channel for ACKs. If the 
source has not yet received an ACK from all M destinations, the packet is retransmitted. This policy 
of relentless retransmissions is assumed throughout the present work. We note that this policy is sub- 
optimal in terms of stable throughput. For instance, in the case of a single destination M = 1, collision 
resolution algorithms such as the one in [9] have been shown to provide a higher stable throughput in 
the infinite-user case. We choose to focus our attention on random access with retransmissions as a first, 
non-trivial step in investigating the stability of random access broadcast. 
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A. Queueing stability 

Let Qn(k) denote the length of the queue at the n th source node at the beginning of the k th slot in 
system S. The evolution of the queue is expressed as follows. 

Qn(k + 1) = (Q„(fe) - B n {k)) + + A n (k), (1) 

where x + = x if x > 0, and otherwise. In the above equation, A n (k) denotes the arrivals to source n, 
where E[A n {k)\ = A n ,Wc and B n (k) denotes completed services (or departures) from n. Thus B n (k) 
takes value 1 if a packet from source n completes service in slot k. We introduce the service rate 
/i„ = linife^oo Pr{B n (k)} as the probability that a packet completes service in the steady-state. 

The vector of queue lengths forms an iV-dimensional, irreducible, aperiodic Markov chain Q(k) = 
(Qi(k),Q 2 (k),. . . ,Q N (k)). The system is stable if, for x € 

lim lim Pr{Q(fc) < x} = 1. (2) 

For our Markov chain Q(k), stability is equivalent to positive recurrence of the Markov chain. We define 
the stability region of the system is the set of all arrival rates (Ai, A2, . . . , A^) for which there exists 
a set of transmission probabilities (pi,P2, ■ ■ ■ ,Pn) such that the system is stable. A primary tool used 
in our work is Loynes' result [10], which tells us that if A n (k) and B n {k) are nonnegative, finite, and 
strictly stationary, then source n is stable if and only if A n < fi n . 

B. Dominant systems 

In the original system S, we cannot easily write down the average service rate \x n of a source because 
the service rate varies depending on whether the other sources are empty or backlogged. Instead, we 
introduce a dominant system <SM which behaves exactly like system S except that all N sources continue 
to transmit "dummy" packets when empty. The dummy packets do not affect the information-carrying 
ability of the source, but their transmission results in a decoupling of the queues. In the dominant system 
e>M, all sources behave as if they are backlogged, the probability of interference from other sources is 
known according to the p n values, and we can easily write down the service rates ■ Let Q^} denote 
the length of the queue at source n in system <SW. It can be shown [2] that if 1$ < fi n then Vx G N, 

Pr{QW >x}> Pr{Q n > x}. 

In other words, the length of the queue in <SM is never shorter than in S. So if we find the conditions 
for stability in then stability in S is implied. Thus, stability in the dominant system is a sufficient 
condition for stability in the original system. 
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C. Throughput Region 

We define the throughput region of the system as the stability region under the assumption that all 
sources are always backlogged. As such, we ignore the burstiness of the arrival process. The throughput 
region of the system S is equivalent to the stability region of the dominant system <SW. As we noted 
above, stability in the dominant system <SM implies stability in the original system S, so the throughput 
region of the system provides an inner bound to the stability region. It has been conjectured (see [11] 
and references therein) that the stability region and the throughput region coincide for N sources and a 
single destination. This conjecture suggests that the zero state Q n = plays no role in determining the 
ergodicity of the Markov chain Q(k). A proof of this conjecture has not been found, and is hindered by 
the fact that the stability region for N sources has not been found. 

III. A NETWORK OF N-2 SOURCES AND M-2 DESTINATIONS 

A. The channel model and service rates 

We introduce a channel model similar to the one in [8] to represent the probabilistic reception and 
MPR that can be attributed to a wireless channel. We define the reception probability when a single 
source transmits as follows. For n,m=l,2, 

Qn\l = Pr{P ac ket from source n is received at destination m | only source n transmits}. 

Additionally, for the N-2 scenario, we will include the possibility of capture or MPR with the following 
reception probability. 

QnXiz = P f {P ac ket from source n is received at destination m | both sources transmit}. 

If both s\ and s 2 are backlogged, the probability that a packet transmitted from s\ is successfully received 
at both destinations is given by 

— (1) (2) , (1) (2) m 

where p^ = 1 — p 2 . In general, r ra ,n=l,2 is the probability that a packet from s n is received at both 
destinations when that source attempts transmission. Similarly, we define 4> n as the probability that a 
successful reception occurs at d^ given that s n transmits and a n as the probability that d^ successfully 
receives when s n transmits. When both sources are backlogged, 4>i and o\ are given by 

0"! = P2^ 1 ) +P2^ 1 ) j2 - (5) 
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We find the service rates and ^2b for the system in which both sources are backlogged by taking 
the expected value of B n . We first condition B n on the receiver state, where r n = {r n l \ r n 2 ^) is a vector 
of binary values indicating whether the packet currently being transmitted by s n has been received at 
each destination. The possible receiver states are r n = {(0,0), (1,0), (0, 1)}. We do not allow for the 
state (1, 1) since upon reaching that state, the source will immediately begin serving the next packet. We 
can write the conditional distribution of B n as given below. 



Pr{B n = l|r n } = ^ 



PnT n , r n =(0,0) 

Pn0n, r n =(0,l) (6) 
Pn&n, r n =(l,0). 

To find the expected value of B n we must first find the steady-state probability of each receiver state. The 
set of receiver states can be modeled by the Markov chain with transition probabilities depicted in Fig. 
12 Let 7To 5 o, tto,i> an d 7Ti o denote the steady-state probabilities of r n when both sources are backlogged. 
Through use of the balance equations for the Markov chain and the equation 710,0 + ^0,1 + ^1,0 = 1 we 
find the steady-state probabilities to be 



7T0,0 



7T0,1 



7T1.0 



+ cr n ){4> n + a n - T n 



+ cr n ){4> n + a n - T n 

4>n(<f>n ~ T n ) 



((f> n + a n )(4>n + &n — Tn) ~ 4> 

Finally, the backlogged service rates // n &, n=l, 2 can be obtained as \i n ^ = E[B n ] = p n r n 7To o + 
Pn4>n^o,i + Pn&nKifi. After simplification, the backlogged service rates can be expressed as 

Pn^n^n((f>n + O n - T n ) 



fJ-nb 



(7) 



We let ni e denote the service rate of s% when S2 is empty and similarly, /X2 e is the service rate of S2 
when s\ is empty. These service rates can be found directly from the backlogged service rates as 



Mle = Ml6 



P2=0 



H2e = M2£> 



Px=0 



(8) 



B. Stability and throughput regions 

We apply the approach introduced in [3] to find the stability region for N-2 and M=2. For fixed 
(pi>P2) the following theorem provides the condition for stability. 
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Fig. 2. The receiver state Markov chain and transition probabilities for N — 2 and M — 2. 



Theorem 1: For a network with N = 2 sources, M = 2 destinations, fixed (pi,P2) and < /u ne , n = 
1, 2, a necessary and sufficient condition for stability is that the arrival rates lie within the union of the 
following two regions. 

A 2 / A 2 \ 

Ai < — Hib +1 /Ule, A 2 < fJ, 2 b- (9) 

^2b V /W 

Al < Ml6, A 2 < — /U 2 fe + ( 1 - — J M2e- (10) 

This theorem is a generalization of the result in [3] and the proof follows the one provided in that work. 
To obtain the stability region over all (^1,^2) we formulate a constrained optimization problem in which 
we fix Ai and maximize A 2 over (pi,p 2 ) subject to Eqns. l9l and [T0l T81. 

Analyzing the throughput region of the system is equivalent to examining a single dominant system, 
and the corresponding condition for stability over all (pi,p 2 ). The stability condition for fixed 
ipiiP2) iS expressed as 

Ai < fx lb , A 2 < fJw (11) 

We can find points on the boundary of the throughput region by again fixing the value of Ai and 
maximizing A 2 over (pi,P2) subject to Eqn. [TT] The following theorem states that the throughput region 
coincides with the stability region; a proof is found in the appendix. 

Theorem 2: For a random access broadcast system with N = 2 source nodes, M = 2 destination 
nodes, and service rates as given in |7]) and ([8]), the throughput region is identical to the stability region. 
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Fig. 3. Stability and throughput regions for N=2 sources and M=2 destinations and two different channel models. The reception 
probabilities for channels (I) and (II) are shown in Table U 



C. Numerical results 

In Fig. [3] we show the stability and throughput regions computed for two different channels with 
reception probabilities given in Table J] The figure demonstrates that as the probability of capture or MPR 
increases, the stability region transforms from a strictly concave region to a convex region bounded by 
straight lines, which agrees with the results in [8] for unicast transmission. These results also demonstrate 
that for a network with N-2 sources and M-2 destinations, the broadcast stability and throughput regions 
coincide, confirming Theorem |2 

IV. A NETWORK OF N SOURCES AND M DESTINATIONS 

A. The channel model and service rates 

The channel model is now a simplified version of the model presented for the N-2, M-2 scenario. 
We assume that whenever two or more sources transmit simultaneously, none of the transmissions are 
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successful. Additionally, we assume that the channel reception probabilities from a source are the same 
for all destinations, q ) = q, = . . . = a, = q„\„, n = 1 . . . N. We refer to the destinations as being 
indistinguishable in this channel model. 

We define the receiver-state variable r n for each source n = 1,2, ... ,N as the number of destinations 
that have received the packet that source n is attempting to transmit, r n E [0,1, . . . , M — 1]. We do 
not allow r n to take value M since as soon as all destinations have received the packet, the source will 
instantaneously revert back to r n = and either begin serving the next packet in the queue of source n 
or become idle if the source is empty. We define the set B as the set of all sources that are backlogged 
at the time source n is attempting to transmit the packet at the front of its queue. Let j3 n denote that 
probability that source n accesses the channel without interference, i.e., 

Pn=Pn[ J] Pi)' < 12 ) 
\leB\n J 

The service process conditioned on r n can be described by 

Pr{£? n = l\r n = r} = (3 n q n \n M ~ r - (13) 

We again develop a Markov chain model for r n as shown in Fig. |4] In this model, transitions "upward" 
can occur between all pairs of states, however, transitions "downward" can only occur between a state 
and the state. Additionally, each state has self-transitions. Let V denote the transition matrix for this 
Markov chain. When source n accesses the channel without collision, which happens with probability f3 n , 
the transition probability matrix will be a matrix V* which depends only on the reception probabilities 
q n \ n . Otherwise, a self-transition occurs, corresponding to a transition probability matrix equal to the 
identity matrix 1. Thus, we can describe V as the convex combination of two probability matrices, 

V = (5 n V* + (1 - (3 n )l. (14) 
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Let 7r be the stationary distribution of V*, iv = ivV*. Clearly iv will also be the stationary distribution 
of V since tvV = f3 n iv + (1 — /3 n )7r. We can solve for tv as follows. Let p* • denote the probability of 
transition in r n from i to j conditioned on s n accessing the channel without collision. The transition 
probabilities are given as 

' (1 " <ln\n) M + < B , ^ = 

< i < M. 



Plo 



(15) 



M- 



Pi,j 



(l-q n \n) M -\ 



(j-i)ln\n^ 1n\n 



\M-j 



i = j 



% < j 



i > j. 



(16) 



In order to satisfy -k; l = V • iTjp* • we have 



i-l 



7T; = ^ TTjP* j + TTjP* i 1 = 1,2,..., M- 1 
j=0 



Ei— 1 
7=0 ^ 



1 - »*• 



(17) 



Together with tt; = 1, we can find the steady-state probabilities tv, which is the stationary 

distribution of V. 

Once the steady-state probabilities of the receiver Markov chain are found, the average service rate is 
given as follows. 

M-l 

Vn = PnJ2 ^n\n M ~ l ( 18 ) 



We further define 



i=0 
M-l 



a, 



i=0 



M-i 



Thus the service rate can be represented in simplified form as 



(19) 



(20) 



This equation summarizes the natural relation between our broadcast problem and the unicast collision 
channel problem [5]. The probability that source n completes transmission of a packet, given by fi n , is 
equal to (3 n , the probability that the source access the channel without collision, times a n , which is the 
probability that all destinations receive the packet conditioned on collision-free access to the channel. In 
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the unicast collision channel problem, we have q n \ n = 1 and a n = 1. Thus for a n = 1 we would expect 
the stability region for our broadcast problem to coincide with the results in [5] for the unicast collision 
channel. This is indeed the case as shown in the following section. 

We define the empty and backlogged service rates as follows. When B contains all N sources, the 
broadcast service rate will take its minimum value fi n b where 



The maximum value of \x n is attained when only source n is backlogged and all other sources are empty. 
We denote this service rate by [i ne 



In deriving bounds on the stability region, we will take advantage of the form of \i n as expressed in Eqn. 
[201 and the fact that /^ n & < \i n < fine- 
Before continuing, we observe the effect of the number of destinations on the stability and throughput 
regions as shown in Fig. [5] These results are for N=2 sources and M = 2,5, and 15 destinations. 
The results are generated using the approach outlined for N-2, M-2 with the exception that we use 
the backlogged and empty service rates for M destinations as given in Eqns. [2T] and |22l The results 
demonstrate that as the number of destinations increases, the broadcast stability and throughput regions 
diminish in size. Additionally, we observe that the stability and throughput regions coincide. This result 
is not surprising since the proof of Theorem [2] holds for arbitrary M with the channel model described 
above. 

B. Stability and throughput regions 

We follow the methodology outlined in [5] to develop bounds on the broadcast stability region. These 
results are a generalization of the bounds on the unicast stability region given in [5]. To begin, by the 
dominant systems argument and Loynes' result, we can develop loose inner and outer bounds on the 
stability region. First, if A n < then since < \i n , queue n must be stable. Furthermore, if \ n < fj, n f, 
for all n, then the entire system is stable. Likewise, if X n > /i n & for all n, then the system is unstable. This 
follows since A ra > /x n & corresponds to instability of all the queues in dominant system «SW, in which 
case all of the queues grow to infinity and the dominant system <SM becomes indistinguishable from the 
original system S [3]. In order to improve upon these loose bounds, we make use of the stability rank 
of the queues as introduced in [5]. Let denote a dominant system in which sources k, k + 1, . . . , N 




(21) 



(22) 
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Fig. 5. Stability and throughput regions for N=2 sources and (I) 2, (II) 5, and (III) 15 destinations. The reception probabilities 
are: = 0.7, g 2 |2 = 0.8. 



transmit dummy packets when empty while sources 1,2, . . . ,k — ldo not. The proof for the following 
theorem is the same as in the original [5] with the exception of the constant a n . 
Theorem 3: Given A at and pn, we order the indices of the sources so that 

M 1 -Pi) < < x n(1 -Pn) (23) 

Then in system S and any dominant system S^ k \ if queue i is stable and j < i, then queue j is also 
stable. 

We can find conditions for the stability of the system through the procedure outlined in [5]. After 
ordering the sources according to stability rank, we first check for stability of queue 1 in system If 
we find that queue 1 is unstable, we can conclude that the entire system is unstable. However, if queue 
1 is stable, we proceed by examining queue 2 in system S^. Queue 1 is known to be stable in due 
to the stability rank. Given stability of queue 1, we will check whether queue 2 is stable in S^. The 
procedure continues in which the stability of queue k in system <S^ is verified assuming that queues 
1, 2, . . . , k — 1 are all stable in S^. If we can finally conclude that queue N is stable in system S^ N \ 
then this implies that the original system is stable. 

Assuming that sources 1, 2, . . . , k — 1 are all stable, we will develop bounds on the average service 
rate ^ for source k in system in order to help determine whether is stable. We can express /J^ 
as 

N 

^ ] =a kPk PP H (1-pj), (24) 
j=k+l 
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\k] 

where P E is the probability that none of the sources 1,2, . . . , k — 1 transmit in the dominant system 
<S[ fe l. One way to bound /J? is by bounding P^' when expressed as 

' = 1 — Pr {only one of sources 1, 2, . . . , k — 1 transmits} 

— Pr {more than one of sources 1, 2, . . . , k — 1 transmit}. (25) 

The bounds on n k result in the following two theorems. 

Theorem 4: Sufficient condition. Given an N source, M destination random access system with \n, 
Pn and the sources ordered according to the stability rank as in Eqn. [23] if Vfc, 1 < k < N, X k < B k , 
where B k is defined below, then the system is stable. 

B k = max(C k ,D k ) 

N 

Bi = a lPl Y[(l- Pj ) (26) 



Cu 



i=2 

akPk 
1-Pk 



na-«)-^-|E(^n(i-»)-^ 

j=k KKk—1 i =1 \ l=fc 



* - r^na-*)(i + g(i-£) » 



i=i \ i=i 



^ y i- Pi J 

(27) 

Theorem 5: Necessary condition. Given an iV source, M destination random access system with 
Pn and the sources ordered according to the stability rank as in Eqn. [23] a necessary condition for stability 
of the system is that Vfc, 

X k <^(U(l-p 3 )-^^-). (28) 

The proofs for Theorems [4] and [5] are described in the Appendix. We make use of these two theorems 
to develop bounds on the stability region by fixing the values of Ai, A2, • • • , Aj\r-i and numerically 
optimizing Aat over all pn € [0, 1]^ subject to the inequalities given in the Theorems. 

The throughput region of the network of arbitrary size can be determined exactly using the backlogged 
service rates in Eqn. [2T] Again, we fix Ai, A2, • • • , Ajv-i and maximize Xn over all pn subject to 

A TV < I^Nb- 
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TABLE II 

Optimization results of stability bounds and throughput for a network of N=4 sources and M=8 and 10 

DESTINATIONS. THE VALUES OF Ai, A2, AND A3 WERE FIXED AND A4 HAS BEEN OPTIMIZED OVER ALL p. 



M 


?i|i 


?2|2 


9313 


?4|4 


Ai 


A 2 


A3 


Stability-upper 


Stability-lower 


Throughput 


8 


0.9 


0.8 


0.7 


0.9 


0.01 


0.01 


0.01 


0.3648 


0.3213 


0.3213 












0.07 


0.02 


0.01 


0.2125 


0.1672 


0.1672 












0.05 


0.05 


0.05 


0.1363 


0.0566 


0.0566 












0.07 


0.05 


0.05 


0.1090 


0.0376 


0.0376 


8 


0.8 


0.8 


0.8 


0.8 


0.01 


0.01 


0.01 


0.2527 


0.2433 


0.2434 












0.07 


0.02 


0.01 


0.1294 


0.1090 


0.1090 












0.05 


0.05 


0.05 


0.0784 


0.0428 


0.0428 












0.07 


0.05 


0.05 


0.0587 


0.0254 


0.0254 


10 


0.8 


0.8 


0.8 


0.8 


0.01 


0.01 


0.01 


0.2329 


0.2236 


0.2236 












0.07 


0.02 


0.01 


0.1153 


0.0951 


0.0951 












0.05 


0.05 


0.05 


0.0651 


0.0318 


0.0321 












0.065 


0.05 


0.05 


0.0503 


0.0196 


0.0196 



C. Numerical results 

A collection of results on stability and throughput for various values of N and M are shown in Tables 
ITTI and |III1 The values in Table JI] demonstrate the effect of the channel reception probabilities and the 
number of destinations on the stability and throughput regions. The trends in these results are the same as 
those observed for a network of N=2 sources. In Table [III] we observe that as the number of sources N 
increases, the stability and throughput regions diminish in size. Furthermore, in all cases, the throughput 
values fall within the upper and lower bounds on the stability values. As such, these results support the 
conjecture that the stability and throughput regions coincide. Of special note, the lower bound for stability 
and the throughput value appear to be equal in many cases. This is not entirely true. In the results shown 
here, the throughput value is in fact slightly larger than the lower bound on the stability value, but the 
difference is at most 10 -6 . 
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TABLE III 

Optimization results of stability bounds and throughput for a network of N=5 and 10 sources and 
M=10 destinations. The reception probabilities are q n \„ = 0.8 for n = l,...,N. The values of Ai, . . . , Ajv_i 

WERE FIXED AND Ajv has been optimized over all p. 



N 


Ai, . . . , Ajv_i 


Stability-upper 


Stability-lower 


Throughput 


5 


[0.010, 0.010,0.010, 0.010] 


0.2078 


0.1939 


0.1939 




[0.070, 0.020,0.010, 0.010] 


0.1051 


0.0789 


0.0789 




[0.035,0.035,0.035,0.035] 


0.0751 


0.0362 


0.0362 




[0.050, 0.035, 0.035, 0.035] 


0.0602 


0.0223 


0.0223 


10 


[0.010,0.010, ...0.010] 


0.1266 


0.0912 


0.0912 




[0.070,0.010, ...0.010] 


0.0679 


0.0252 


0.0252 




[0.017,0.017, ...0.017] 


0.0621 


0.0137 


0.0137 




[0.020,0.017, ...0.017] 


0.0591 


0.0108 


0.0108 



V. Conclusions 

In this work we have investigated the information-carrying ability of source nodes which compete 
for access to a shared channel in order to broadcast messages to a common set of destination nodes. 
The scenario we consider is relevant to broadcast applications and network discovery and control. We 
formulated the problem as a network of queues which interact through the shared channel and investigated 
the stability and throughput regions of the source nodes. Our results strengthen an unproven conjecture 
that the stability and throughput regions of finite-user random access coincide. We also demonstrate that 
the throughput and stable throughput regions diminish as the number of source or destination nodes 
grows. Our intention in this work is to present results relating to the scaling laws for wireless networks, 
which often show a degradation in throughput as the number of nodes increases. Our results show a 
similar trend from a multicast stability point of view. 
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VI. Appendix 

A. Proof of Theorem [2] 

We show that the boundaries of the stability and throughput regions, given by the results of two 
constrained optimization problems, are identical. We replace Ai by x and A2 by y. Furthermore, we 
express the backlogged service rates in Eqn. |7] as follows. 

Mib =Pi/i(P2), M26 =V2h(vi)- (29) 

The empty service rates can then be expressed as /xi e = pi/i(0) and /i 2e = ^2/2(0). In these expressions, 
/i(p 2 ) is functionally independent of p\ and decreasing in p 2 , corresponding to fi\b < \i\ e f\{p2) < 

/i(0). Similarly, / 2 (pi) is independent of P2 and decreasing in p\. With this notation in place, the 
boundary of the region given in Theorem Q] for fixed (pi,p2) can be written as follows. 

x = pi/i(0) — , 0<y<P2f2(Pi) (30) 

P2MP1) 

f , n s Xp 2 (/ 2 (0) - /2(pi)) n . . , , x 

2/ = P2M0) tV^; , < X < P1/1CP2) (31) 

Plfl{P2) 

To find the stability region, we should maximize the expressions in (l30l and (OTT ) over (^1,^2) an d 
take the intersection of the regions bounded by the resulting curves. We note that this is not a standard 
optimization problem because the objective function is piece-wise and non-differentiable at a point in its 
domain. An example of an analytical solution to this optimization problem is given in [8] for the single 
destination case. 

The boundary of the throughput region for fixed (pi,P2) as given in Eqn. [TT] can be written as follows. 

x = Pifiipz), < y < p 2 / 2 (Pl) (32) 
V = £2/2(^1), < x < pifi(p 2 ) (33) 



The throughput region is found by maximizing the expressions in (132] ) and (1331 ) and taking the intersection 
of the resulting regions. Consider Eqn. ( T33T ) in which we wish to maximize y over (pi,p 2 ). Note that the 
constraint x < Pifiipz) is a lower bound on p\ over which we perform the maximization. Since / 2 (pi) 
is decreasing in pi, the lower bound p\ > x//i(p 2 ) provides an upper bound on y, and this upper bound 
can be achieved when maximizing (|33l ) over (pi,p 2 ). Then maximization of y in ( f33T > is equivalent to 
maximizing y as follows. 



V = P2f2(pi) 



, 0<x<Pifi{p2) (34) 

Pi= 7T^y 
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To see that this is identical to (I3TT) . we write /2Q31) as follows. 



/ 2 (pi) = / 2 (0)-(/ 2 (0)-/ 2 (pi)) (35) 

= / 2 (0)- Pl(/2(0) " /2(Pl)) (36) 
Pi 



Then 

and the maximum y in ( f33T > is identical to the maximum y in (|3T| ). Similarly, the maximum of x in (|32 
is identical to the maximum of x in (l30l). 



= P2/2 (0)-^iM^M (37 ) 



£. Proof of Theorem |4] 

Our proof follows the one given in [5]. A sufficient condition for stability corresponds to bounding 

\k] 

fi k from below. We obtain two separate lower bounds, Ct and D^, and take their maximum to obtain 
the result. The bound Ct is derived through Eqns. |24] and [25] Since sources 1, 2, . . . , k — 1 are stable, 
the following holds. 

k-l 



Prjsuccess by one of the sources 1, 2, . . . > k — 1} = Aj 

i=i 

> ( min ai ] Pr {only one of the sources 1, 2, . . . , k — 1 transmits} JT(1 — Pj) (38) 
This indicates that 

V fc_1 A- 

Prjonly one of the sources 1, 2,..., k - 1 transmits} < -. \ i=1 1 . (39) 



For the next term in P l E J we have 



Pr{a packet from j collides with others from 1, 2, . . . , k — 1} 

= -Pr{j transmits} — Pr{j transmits & no others from 1, 2, . . . , k — 1 transmit} 

Pr{j transmits & no others transmit} 



PjPrjj backlogged} 

lli=fc(l - Pi) 

_Aj A 



Pi^Sr-— ^ (4°) 
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h'l [fel \k] 

We also make use of the result shown in [3] that for j < k, fij < fij . Thus, the second term in P E 
is upper bounded as follows. 

fc-l 



Pr {more than one of the sources 1, 2, . . . , k — 1 transmits} < - ^ ( pj 

2 i=A ^ 

By combining Eqns. |25l [39l and @j] we obtain the lower bound 



-\j ^3 ] 

(41) 



^^-pi^y*^-^' (42) 



Together with Eqn. |24l this provides Cfc, our first lower bound on /il 

[fel 

The other lower bound on fj,i is derived using an approach from [3]. We have adapted it for the 
broadcast problem below and refer to it as D^. 

[fel a kPk 

u, > PWno other source transmits [ 

k l-Pk 

H Pr{one of the first 1, 2, . . . , k — 1 sources transmits but is empty} 

1 ~Pk 

Note that we can express the exact value of ^ as ^ = aiPin^=2(l — Pj)- beginning with 
B\ = ^ we can iterate through k values 2, . . . , N to obtain the result. 

The proof of Theorem [5] develops upper bounds on /^Jj. an d thus on Pjp . The technique is similar to 
the one outlined above in finding Ck- 
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